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Abstract. The maximal group of Lie point symmetries of a system of 
nonlinear equations used in geophysical fluid dynamics is presented. The 
Lie algebra of this group is infinite-dimensional and involves three arbitrary 
functions of time. The invariant solution under the rotation and dilation is 
constructed. Qualitative analysis of the invariant solution is provided and 
the energy of this solution is presented. 
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1 Introduction 



This is a continuation of the paper [1]. We present here the Lie algebra of 
the maximal group of Lie point symmetries for system nonlinear equations 

A-0t - gp x - fv z = ip x Aip z - 4>zMx , (1.1) 

Vt + flpz = i>xVz ~ IpzVx , (1-2) 

N 2 

Pt + —i>x = i>xpz - *Pzpx (1.3) 

used in geophysical fluid dynamics, e.g. for investigating internal waves in 
uniformly stratified incompressible fluids (oceans). Here g,f,N are con- 
stants and A is the two-dimensional Laplacian: 

A = Dl + Dl- 



2 Symmetries 
2.1 General case 

The point symmetries of Eqs. (1.1)-(1.3) have been computed with the help 
of DIMSYM 2.3 package. The maximal admitted Lie point transformation 
group is infinite for arbitrary constants / and N. If / ^ 0, the group is 
generated by the infinite-dimensional Lie algebra spanned by the following 
operators: 



Xi = w 



d_ 

dv 



X, = 



X 7 = 



b(t) 
c(t) 



x 2 


d 

dp ' 


' d 


-/-I 

dv 


dx 


r 9 


n 2 d 



d 



at' 



dz ^ g dp 



- c'(t)x 



d_ 

dip 



(2.1) 



Xo = 



d d d d n , d 

+ Z 1T + V 7T + PIT + 2 ^^T ' 
dx dz dv dp dip 

d d d , d „ d N 2 d 
t— + 2x— + 2z— + 3tp— - 2 fx— + 2—z—, 
dt dx dz dtp dv g dp 

4-4-7 {9P+(f 2 -^}l + l .ifv+(f 



dx dz f " dv g 

Here a(t), b(t) and c(t) are arbitrary functions of time t 
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Remark 2.1. The presence of the arbitrary functions a(t), b(t),c(t) in the 
symmetry Lie algebra is a characteristic property of incompressible fluids 
([2], see also [3]). Namely, the operator X3 generates the group transforma- 
tion -0 = ip + £ 3 a (t) °f t ne stream function -0, where £3 is the group parame- 
ter. The invariance of fluid flows under this transformation is quite obvious 
because the velocity vector (ip z ,v, —ifj x ) is invariant under this transforma- 
tion. The operators X 5 ,X 6 express the invariance under the generalization 
x = x + E5b(t), z = z + Sq c(t) of the coordinate translations and the 
Galilean transformations. They provide a generalized relativity principle 
for the Euler equations in terms of conservation laws (see [4], Section 25.3). 



2.2 The case f = 

In order to include the special case / = 0, we multiply the operator X 9 by 
the constant / and consider the operator 



K = f 







o 



x- 



dx dz. 

Then we let / = and obtain the operator 

X>=-[gp-N 2 z] 



[9P+(f 2 -^)z}^ + l[fv + (f-N 2 )x}§- p 



d_ 

dv 



admitted by Eqs. (1.1)-(1.3) with / = 0. The solution of the determining 
equations shows that X' 9 is a particular case of a more general symmetry 
involving an arbitrary function of two variables. Namely, the system (1.1)- 
(1.3) with / = admits the infinite-dimensional Lie algebra spanned by the 
following operators: 



d 



d 



X 1 = h(v,gp-N i z) ir , X 2 = —, X. 



dv 



dp 



X 6 = c(t) 



d_ 

dz 



d d 
dx dz 



n 2 d 

g dp 
d 



// \ d 



(2.2) 



d 



P- 



dv dp 



2if; 



d_ 

dip ' 



9 n d n d n , d „N 2 d 
X 8 = t— + 2x— + 2z— + 3ifj— + 2—z—, 
dt dx dz dip g dp 

where h(v, gp — N 2 z) is an arbitrary function of two variables. The operator 
X\ in (2.1) is obtained from the operator X 1 in (2.2) by taking h — 1. 



3 



3 Invariant solution based on rotations and 
dilations 



3.1 The invariants 

We will investigate here the invariant solutions with respect to the dilations 
and rotations with the generators X 7 and X 9 . Let us introduce the notation 

v* = fv, u = gp, a = f 2 -N 2 (3.1) 

and write the operators X 7 , X 9 in the form 



(3.2) 



v d d d d d 

X l = X TT + Z TT + U TT + V *T> h 2 ^^T ' 

ox oz ou ov* oip 

d d d d 
X 9 = z- x— + k + ai)- (u + az) —— ■ 

ox oz v ou ' ov* 

The operators (3.2) coincide with the operators (3.17) from [5], 

d d d d.d 
ox oy ou ov ow 

d d d d 

X 2 = x— + (v + ax + (3y)- (u - fix + ay)— , 

ox oy ou ov 

with k = 2 and (3 = upon identifying v with v* and y with z. Hence, a basis 
of invariants for the operators (3.2) contains the time t and the invariants 
(3.20) from [5] which have now the form 

1 

Ji = 77 (xu + zv* + axz) , 

X 2 + z l 

J 2 = 2 1 2 ( XV * ~ ZU + 77 ~ z2 )) ' 

x A + u V 2 / 



J. 



x A + y z 
x 2 + y 2 



It is more convenient for our purposes to use, instead of these invariants, 
the equivalent equations (3.19) from [5] which are written now as follows: 

u = J x x - ( J 2 + ^ z , 

u* = Ji z + ( J 2 - |) a , ( 3 - 3 ) 
ij = {x 2 + z 2 ) J 3 . 
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3.2 Candidates for the invariant solution 



Knowledge of a symmetry algebra allows one to obtain particular exact so- 
lutions to differential equations in question. These kind of solutions were 
considered by S. Lie [6]. They are known today as group invariant solu- 
tions (briefly invariant solutions) and widely used in the modern literature, 
particularly in investigating nonlinear differential equations. 

The general form of regular invariant solutions is obtained from Eqs. 
(3.3) by setting 

J! = R(t), J 2 = V(t), ,h = 4>{t) 

with undetermined functions R(t), V(t), <f>(t). Invoking the notation (3.1) 
we arrive at the following general form of candidates for the invariant solu- 
tion with respect to the dilations and rotations with the generators X 7 and 
X 9 from (2.1): 



/ 
1 

P = ~ 
9 



N 2 — f 2 ' 
R(t) z + V(t) x + x 

N 2 - f 2 " 
R(t) x - V(t) z + — — z 



(3.4) 



i{; = (x 2 + z 2 ) ( J ) (t). 



Remark 3.1. Solving the Lie equations for the operator X 9 from (3.2) and 
using the notation (3.1), one can verify that the operator X 9 from (2.1), 



X n = z^--x^--\ [g P + (f - N 2 )z] ^ + ~ g [fv+ (f - N 2 )x] 9 



dx 



dz f 



dp' 



generates the following one-parameter transformation group with the pa- 
rameter e : 



x = x cos e + z sin e, z = z cos e — x sin e, 
gp = gp cos e + fv sin e — (N 2 — f 2 ) x sin e, 
fv = fv cos £ — gp sin e + (N 2 — f 2 ) z sin e, 
i = t, ijj = ip. 



(3.5) 



One can verify by inspection that the transformations (3.5) leave invariant 
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Eqs. (3.4): 



1 

s = 7 



N 2 - f 2 ' 

R(t) z + V(t) x + — x 

2 

N 2 - f 2 " 
R(t) x - V(t) z + — — z 



ij = (x 2 + z 2 ) ( f ) (t). 



3.3 Construction of the invariant solution 

It remains to determine the functions R(t), V(t), <f>(t) by substituting the 
expressions (3.4) for p, v, if) in Eqs. (1.1)-(1.3). 
Differentiating (3.4) we obtain: 



v t = j[R'z + V'x}, v x = j 



AT 2 -/ 2 



+ V 



v, — —R, 



Pt — -[R'x — V z] , p x — -R, p z = - 
9 9 9 



N 2 -f 2 



f 



V 



(3.6) 



rjH = (x 2 + z 2 ) (j>\ rj} x = 2x (f>, ij> z = 2z (f>, A^ t = 40'. 
Substitution of (3.6) in Eqs. (1.1)-(1.3) yields: 
2<j/-R = 0, 

[V - 2R<P] x+[R' + 2V(f) + (N 2 + f 2 )<f)} z = 0, 
[R' + 2V(j> + (N 2 + f 2 )(f)} x-[V- 2R(j)\ z = 0. 

Since V, R, (p depend only on t, Eq. (3.8) implies that 

V - 2R4> = 

and 

R! + 2V(f) + (N 2 + f 2 )(f) = 0. 



(3.7) 
(3.8) 
(3.9) 

(3.10) 
(3.11) 



Eq. (3.9) is satisfied due to Eqs. (3.10), (3.11). Hence, Eqs. (1.1)-(1.3) are 
reduced to Eqs. (3.7), (3.10), (3.11). 
Let us write Eq. (3.7) in the form 



R = 2(j)'. 



(3.12) 
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Substitution of the expression for R into Eq. (3.10) yields V = ', 
whence upon integration 



V = 2(f) 2 + A, A = const. 



(3.13) 



Finally, substituting Eqs. (3.12) and (3.13) in Eq. (3.11) we obtain the 
following nonlinear second-order ordinary differential equation for <fi(t) : 



4>" + 20 3 



A 



f 2 + N' 



= 0. 



(3.14) 



Thus, we have arrived at the following result. 



Theorem 3.1. The solutions of the system (1.1)-(1.3) that are invariant 
with respect to the dilations and rotations with the generators Xj and X 9 
from (2.1) are given by 



1 

v = l 



24> 2 {t) + A + 



N 2 -f 



2<f>'(t)x - 2(f) 2 {t) + A- 



x + 2<f)'(t)z 



N 2 -f 



(3.15) 



iP = (x 2 + z 2 )<P(t), 



where 4>{t) is defined by the differential equation (3.14) and A is an arbitrary 
constant. 



3.4 Qualitative analysis of the invariant solution 

One can integrate Eq. (3.14) once, e.g., upon multiplying by 2(f)' and obtain 

/2 + 4 + ( A + f 2 + N2 \ 02 = congt ( 3 16 ) 



We will analyze the behavior of the solutions to Eq. (3.16) under the as- 
sumption that the expression in the parentheses is a non-negative constant 
which we denote by K : 

f 2 + N 2 

K = A+ J —- , K>0, (3.17) 

and write Eq. (3.16) in the form 

4>' 2 + 4 + K(f) 2 = B 2 , B = const., (3.18) 
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or solving for 0' : 



0' = ±^B 2 - 4 - K(j) 2 . (3.19) 

Note that <f>{t) = solves Eq. (3.14). Let us turn to Eq. (3.19). When 
is small, i.e. close to the trivial solution 0(t) = 0, then 

B 2 -<f- K<j? » B 2 

and hence 0' is close to the constant value 

0' « ± B. 

When 4>(t) varies according to Eq. (3.14), then |0'| decreases since 

B 2 - 4 - K0 2 < B 2 
when 7^ 0. We obtain 0' = when 0(t) = C*, where 

^ 9 -K + VB 2 + K 2 

Cl = \ (3.20) 

If |0| > | C* |, then B 2 - 4 - Ktp 2 < 0, and hence Eq. (3.19) does not have 
a solution. We have arrived at the following significant results. 

Theorem 3.2. Provided that the condition (3.17) holds, the solutions of 
Eq. (3.19) are bounded oscillating functions 0(i) satisfying the condition 

-C* < 0(t) < C*, (3.21) 

where C* is the positive constant defined by Eq. (3.20). In this notation, 
the invariant solution (3.15) is written as follows: 

v = j[(2 ( p 2 (t) + K-f 2 )x + 2^(t)z], 

p = I [2<f>'(t)x - (20 2 (t) + K - N 2 ) z] , (3.22) 
rj,= (x 2 + z 2 )<p(t). 

Remark 3.2. The invariance of the solution (3.15) with respect to rota- 
tions (rotational symmetry) means that it has the same values on any circle 

2,2 2 

x + z = r 

with a given radius r. The invariance under dilations means that we can 
obtain the solution at any circle just by stretching the radius r. According 
to Theorem 3.2, this solution is given by bounded oscillating functions. 
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4 Energy of the rotationally symmetric so- 
lution 



The conservation of energy for Eqs. (1.1)-(1.3) has the form [1] 

N 2 



d 
~dt 



v 2 + ^ P 2 + m\ 2 



dxdz = 0. 



(4.1) 



Hence, the energy density is 



E = v 2 + ^p 2 + m\ 2 . 

For the rotationally invariant solution (3.22) we have 

|VVf =4{x 2 + z 2 )<p 2 {t). 



(4.2) 



(4.3) 



Substituting the expression (4.3) and the expressions (3.22) of v and p in 
Eq. (4.2) we obtain the following energy density for the invariant solution 
(3.22): 



e =4 (i - ^) (x 2 - z 2 ) im + K ] m + (f-j) 2 * : 

+ (N - I ) 2 z 2 + 4 (1 - ^) [20 2 (t) + AT] 0'(t) . 



(4.4) 



22 April 2009 
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